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1 Introduction 

The debate about the meaning of the wave function in quantum mechanics is 
still going on, seventy years after the birth of the Copenhagen interpretation Q. 
Despite the enormous successes achieved by quantum mechanics and quantum 
field theory, we are not yet able to answer the following simple question: should 
we consider the wave function as a real objective field, or merely as a useful tool 
for computing probability amplitudes? 

There are (at least) two obstacles to a realistic interpretation of the wave 
function: 

1. For a N-body system, the wave function is defined in a 3N-dimensional 
configuration space, while we would expect a "real" field to be defined in or- 
dinary 3-dimensional space. As a consequence, entanglement between states 
of distant systems is allowed, and this gives rise to the counter-intuitive 
behaviour known as "quantum non-locality" . 

2. The wave function spreads out with time. This property is strictly related to 
the Heisenberg uncertainty principle: if a particle is well localized in position 
space at a certain time, then it will be poorly localized in momentum space, 
and therefore the dispersion of the wave packet will grow with time. Our 
experience with macroscopic objects tells us that they do not spread out with 
time, but in this case the quantum mechanical explanation may be that for a 
macroscopic object the time needed to obtain an observable dispersion of its 
wave packet is greater than the age of the universe. However, microscopic 
objects too are always detected in small space regions and it is hard to 
believe that immediately before their detection they were spread out at long 
distances from the detection point. 

In this paper we conjecture that maybe the two above-mentioned problems 
are a consequence of the particular choice of representation for the wave function: 
with a suitable change of representation, it is possible to find solutions to the 
wave equations which do not spread out with time. We'll see that this result can 
be obtained at the cost of introducing a second space-time coordinate xd) we 
will consider this new coordinate as an "auxiliary" coordinate, in addition to the 



"physical" coordinate xs- Roughly speaking, the wave function will be spread out 
in the direction, while being perfectly localized in the "physical" xs direction. 

We will first examine the non-relativistic Schrodinger equation (Section |2|) and 
then the Dirac equation (Section ||), both for free fields and for interacting fields 
(we will focus on the electromagnetic interaction); finally (Section [|) we will draw 
some conclusions and propose some possible interpretations of the results obtained 
in Sections ^ and ||| 



2 The non-relativistic Schrodinger equation 

Let's consider a quantum system describing a non-relativistic particle moving in 
one dimension under the effect of an external potential V(x). If we choose the 
basis formed by the eigenvectors \x) of the position operator Q, then the repre- 
sentation of an arbitrary state is given by the wave function ip (x) = (x\ip) and 
the Schrodinger equation is 

The time evolution of the system is a unitary transformation, and can be written 
shortly as ip(x,t) = e~n m tp (x, 0). Besides, the position operator Q and the 
momentum operator P are hermitian operators defined by 

Q^j = x%b Pip = -ih^ (2) 

ox 

and they too generate unitary transformations of the form e ta( °^ and e JaP ; the 
transformations generated by the momentum operator are simply the space dis- 
placements, i.e. e l ^- xP / h iff (x) = if> (x + Ax). 

We now define a new field (p (x, y), which depends on the wave function tfj (x) 
through the following relation: 

ip (x, y) = if) (x) if)* (y) (3) 

If we apply an arbitrary unitary transformation to the wave function ip(x), the 
field ip(x,y) will undergo a corresponding transformation. Through this corre- 
spondence, we can obtain the form of the generators in a new representation of 
the group of unitary transformations. In the case of the hamiltonian operator H, 
which generates the time evolution of the system, we easily obtain: 

Equation (Q) is equal to the Schrodinger equation for density matrices, which 
defines the time evolution for mixed states, and therefore we will sometimes refer 
to the new representation as to the "density matrix representation" ; in this paper, 
however, we will try to interpret the field ip as a real objective field, rather than 
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just a tool for computing probabilities. As for the generators Q and P, in the new 
representation we obtain the expressions 



f dip dip \ 

The definitions (|j) and (^) can be used for an arbitrary complex field (p(x,y), 
not only for the fields expressed by the relation (||). However, we will consider 
only fields for which 

ip (y, x) = ip* 0, y) (6) 

This condition is needed to obtain real values for the observable quantities, as 
shown by expressions (|T^)-(|l5|) below. From (|6|) it follows that we could use a real 
representation for the field ip, simply by defining the real field 

(p R = Re{ip} + Im{tp} (7) 

Indeed, Re {ip} is symmetric and Im {</?} is antisymmetric with respect to the 
permutation x <-> y, and therefore they can be unambiguously obtained from 
the real field ipn. However, in the rest of this section we will use the complex 
representation because it produces simpler expressions. 

Now we want to define in the density matrix representation the observable 
quantities, which are the position, the momentum and the energy. In the original 
representation the natural definitions are the following: 



Q = J V* (x) x V 0) dx (8) 

p = J ^* ( x ) (-ibjA V 1 (x) dx (9) 
f ( h 2 d 2 \ 

E = J r (x)i + v ( x ) J ^ ( x ) dx ( 10 ) 

In the probabilistic interpretation of (jl]), the above definitions are simply the mean 
values of the corresponding quantum observables. However, if we consider (|l]) as 
a classical field equation and set V(x) = 0, then it is easily seen that (0) is the 
conserved quantity associated to the space invariance of the system, (|10|) is the 
conserved quantity associated to the time invariance of the system, while (|8|) is 
the center of mass, since @ and (^) satisfy the following classical relation 

Q(t) = — t + Q(0) (11) 
m 

Besides, (|l0|) is exactly the hamiltonian functional needed for applying to ((!]) the 
canonical hamiltonian formalism. Thus the definitions (||), (@) and (10) for posi- 
tion, momentum and energy are reasonable also within a classical interpretation 
of the Schrodinger equation (jlj). 
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For extending the above definitions to our new representation it is convenient 
to introduce a new pair of space coordinates 



x s 



2 ( x + y) 



xd = v 



X 



(12) 



In the rest of this section we will write loosely <p (x, y) or ip (xs,xd), meaning the 
same field expressed in two different coordinate systems. By means of the relation 
(0) we obtain the following expressions in the new representation: 



Q 



x s ip(x s ,x D ) 



x D =0 



dx s 



ih 



d(p(xs,x D ) 



dx 



D 



dx s 



x D =0 



E 



h 2 d 2 ip(x s ,x D ) 



2m 



dx 2 



+ V(x s ) tp(x s ,x D ) 



D 



(13) 



(14) 



(15) 



x D =0 



It can be shown that P and E are still conserved in the case V{x) = 0, while in 
the same case the classical relation ( |TTD is still satisfied. 

Thus the physical quantities depend only on the value of the field cp (and its 
space derivatives) in the region where xd = 0, or x = y. Therefore we are induced 
to consider xs as the "physical" space coordinate, while xjj will be considered 
as an "auxiliary" coordinate, which has physical effects only around the point 
xd = 0. Besides it can be shown that the momentum operator P, as defined in 
(||), satisfies the following relation: 



J Ax P/h 



if (xs, xd) = <p (x s + Ax, x D ) 



(16) 



Since we already know that the momentum operator generates the spatial displace- 



ments, expression (16) adds more evidence to our belief that xs is the "physical" 
position coordinate. 

Now we want to express in the density matrix representation the unitary trans- 
formation which allows to switch from the position space to the momentum space. 
We know that this transformation for the original Schrodinger equation is given 
by the Fourier transform 



1 



ip (x) e~ ikx dx 



(17) 



We will again use the fundamental relation @, exploiting the property that it 
commutes with all unitary transformations, and then extend the result to an 
arbitrary field <p. The final expression is 



$>(k x ,k y ) = — I I tp(x,y)e 



From 



*2/y — 2~ / / w^y^ vV dxdy 
we easily find that in momentum space the operator P is given by 



(18) 



P ^ = (p x - Py) ^ 



(19) 
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where we set p x = hk x and p y = hk y ; we note that expression ( |l9| ) has exactly the 
same form as the definition (|5|) for the position operator Q. Besides, the quantity 
P defined by (|14|) is expressed in momentum space by 



where we have defined 



k s 



P = / p S $(ks,k D ) 



1 



(k x + ky) 



dk s 



kf) — ky k x 



(20) 



(21) 



while obviously ps = hks- Again we note that expression fl20| ) has exactly the 
same form as the definition (|l3|) of the quantity Q in position space. This formal 
equivalence between the operators P and Q, and between the observable quantities 
P and Q, extends to our new representation the simmetry between position space 
and momentum space which is well known in non-relativistic quantum mechanics. 

In the free particle case, where V(x) = 0, the time evolution of the system in 
momentum space is simply given by 



at 



2m 2m 



(22) 



while the energy E becomes 



2 

E= [ ^-<f>(k s ,k D ) dk s 
J 2m k D =0 



(23) 



From ( p0| ) and (23) we are again induced to consider ks as the "physical" momen- 
tum coordinate, while ku will be considered as an "auxiliary" coordinate. 

Now we will answer a very interesting question: what happens to the Heisen- 
berg uncertainty principle in the density matrix representation? In the original 
representation the Heisenberg principle follows from the fact that the operators Q 
and P do not commute: 



[Q,P] = QP-PQ = in (24) 

On the contrary, it can be easily shown that in our new representation the op- 
erators P and Q, as defined by (||), do indeed commute. So it seems that the 
Heisenberg principle is not true in this representation, and indeed we will now 
show that there exist fields ip which are perfectly localized both in position space 
and in momentum space: let's consider the field 

<Po(xs,x D )=5(xs-x )e- ikox ° (25) 

Having identified xs as the "physical" position, we can say that ifo is perfectly 
localized at the position xq, since we have ifo = for xs ^ xq. Applying the 
transformation (]i~8| ) to ipo we obtain in momentum space 

*o(k s ,k D ) = 5(k s -k )e ik ° xo (26) 
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and again, having identified k$ as the "physical" momentum coordinate, we can 
say that $o is perfectly localized in momentum space, with momentum po = hko. 
Besides, applying to <po the expression dl3|), we indeed obtain Q = xq, while 
applying to (fo the expression ( |l4|) we obtain P = fi ko = po as expected; the 
same value for P is obtained in momentum space, i.e. applying the expression 
© to $ - 

In the case of a free particle, where V(x) = 0, we easily find that the 
energy, obtained from fll5| ) or from (|23|), has the correct classical value, i.e. 
E = f^k^jlm = p\j2m. But there is a far more interesting result: by setting 

x = v t k = — = —— (27) 

which are the classical position and momentum for a particle moving with constant 
speed vo, we find that ( p5| ) is a solution of the motion equations (|j). Thus ( |25| ) and 
(p7|) define a field ipo(xg,XD,t) which can be interpreted as a non-relativistic point- 
like particle moving at constant speed vq with position and momentum perfectly 
known at every time t. This result, which seemed incompatible with quantum 
mechanics, has been obtained simply by carrying the Schrodinger equation (||) to 
the density matrix representation. 

Now we turn to the general case of a potential V (x) ^ 0. We are interested 
in finding the stationary states, i.e. those states which are sinusoidal functions 
of the time t. In the original representation, the stationary states ip a were the 
eigenfunctions of the hamiltonian 

ih^=K?p a = E a ip a (28) 

while the set of all eigenvalues E a was the energy spectrum of the system. The 
same is true in our new representation, but now the hamiltonian H is expressed 
by equation (||). It is not difficult to show that in our case the eigenfunctions ip a p 
and the eigenvalues E a g are simply given by 

<^q/3 {x, y) = ip a (x) ipp (y) E aP = E a - Ep (29) 

Thus the natural frequencies in our representation correspond to the difference 
between two energy levels of the original quantum system. But we know very well 
that the individual energy levels are not really observable; the only observable 
quantities are exactly the "jumps" between different energy levels expressed by 
(p9"D- Therefore we can say that we have reproduced the same energy spectrum 
as the original quantum system; besides, we have eliminated the ambiguity about 
the energy zero-point which is typical of all quantum systems. 

In conclusion, simply by changing the representation of the Schrodinger equa- 
tion we have succeeded in bringing together two concepts which are widely held 
as incompatible: on one hand we describe point-like particles, with position and 
momentum perfectly known at every time t, on the other hand we allow the exis- 
tence of discrete energy spectra. Besides, in the case V(x) =0 we can add together 
an arbitrary number of solutions of the form (|25|), thus obtaining a system of non- 
interacting particles; the field describing this system depends always on the same 
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two variables xs and xd, so the dimension of the space on which the field is de- 
fined does not change with the number of the particles. Finally, the transition 
to our new representation involves the definition of an "auxiliary" position coor- 
dinate Xd-, in addition to the "physical" coordinate xs- The existence of extra 
space-time dimensions is obviously not an original idea: it was first proposed by 
Kaluza and then developed by Klein @. In recent years, I find very interesting the 
attempt made by Hasselmann |J to derive all elementary particles properties from 
Einstein's gravitational field equations in a higher-dimensional matter-free space. 
Our model has in common with the cited works the periodicity of the solutions 
with respect to the extra space dimension, as shown by (p5|). 



3 The Dirac equation and the electromagnetic inter- 
action 

We now try to extend the results obtained in Section ||] to the case of a relativistic 
wave equation. We choose the Dirac equation, rather than the Klein-Gordon 
equation, for the simple reason that all known elementary particles are half-spin 
fermions, while no elementary spinless particle is known to exist. We start from 
the non-covariant form of the Dirac equation, which separates the time derivative 
from the space derivatives: 

i^ = (-ia k d k + m(3)7p = K iP (30) 

where the implicit sum is over k = 1,2,3. 

The operator Ho defined by equation ( p0| ) generates a unitary transformation 
on the spinor field vp(x.), where now ip depends upon a three-dimensional space 
coordinate x; we then follow the same method already seen in the previous section 
and define a new representation of the group of unitary transformations by means 
of the fundamental relation: 

y>(x,y) =V(x)V f (y) (31) 
The time evolution in the new representation is then given by 



i — = R p = -i (a k d Xk ip + d yk ipa k ) + m((3p- p (3) (32) 



We now extend equation (32) to an arbitrary 4x4 complex matrix <p (x,y); how- 
ever, as in the previous section, to obtain real values for the observable quantities 
the anti-hermitian part of the matrix field p> must vanish, and therefore we will 
only consider fields pi (x, y) satisfying the relation 

ip (y, x) = p ] (x, y) (33) 

Again we could use a real representation by defining ip r = Re {p} + Im {p}, but 
we prefer the complex representation because it produces simpler expressions. 



We are now able to find solutions to the motion equations (32) which are 



perfectly localized in physical space and do not depend on time; these solutions 
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are thus good candidates to represent point-like particles at rest. To obtain this 
result, we first define a new pair of space coordinates 

1 



x s = ~ ( x + y) 



XD 



X 



and then choose a particular representation for the matrices (3 and a k 



I 
-/ 



o k 
a k 



(34) 



(35) 



where a k are the well known Pauli matrices, while / is the 2x2 identity matrix. 
Our solutions are then given by: 



(36) 







4m6 


1836 i 816 + 826 




1 













¥A = 


4m 


-i 836 














d 2 5-idi8 














" 












1 


Am6 i 


did — 826 


-i 836 


<PB = 


4m 


-idi5-d 2 5 














id 3 6 



















" 






1 





d?5-i&5 




= 


4m 


id^S i 816 + 826 


—4m 6 





























-i 


8 X 6- 


82 




1 








1836 


fD = 


4m 

















1816 — 826 —1836 


— 4mi5 



(37) 



(38) 



(39) 



where we wrote simply 6 for <5(xs). If we identify, as in the previous section, 
the physical position with the coordinate xs, then the solutions (|36|)-(|3"9|) are 
perfectly localized in the physical position x§ = 0. At this point we would like to 
find solutions corresponding to particles moving at constant speed; these solutions 
will be easily obtained from (^)-(|39]) by means of Lorentz transformations, but 
first we must investigate the Lorentz covariance of our new representation. 

There is only one way to recover a formal Lorentz covariance in our new rep- 
resentation: since special relativity treats position and time as components of the 
same vector x^, and since we have duplicated the position coordinate x, we in- 
evitably must duplicate also the time coordinate t. Thus we are led to introduce 
two time coordinates t x ,t y and a second pair ts,tz> obtained through a relation 
analogous to (34). In expressions ([3l]) and (|33| ) then the three- vectors x and y 
must be replaced by the four- vectors x = (t x ,x) and y = (t y , y), while the motion 



(40) 



equations (^2|) are extended to the following covariant form 

„ 8 



8x^ 



<p{x,y) -m<p(x,y) = 



8 



which is formally identical to the covariant Dirac equation in the usual represen- 
tation. Equation (|40|) can be written in the equivalent form 



d 



(41) 



which is obtained by taking the adjoint of (^) and by exchanging x with y. If we 
separate the time derivative from the space derivatives, (^0|) and ( f4l|) become 



. dip 
dt x 

dip 
dt, 



-i a k d Xk ip + mf3 ip 
-id Vk ipa k -m<p(3 



(42) 
(43) 



where we applied the usual definitions /? = 7° and a k = ^/° j k . By adding together 



(42) and (43), we obtain the derivative with respect to ts = (t x + t y ) /2, which is 
equal to the motion equations ( |32| ) in the case t x = t y = t; by subtracting (42) 



from ( f43| ) and dividing by two, we obtain the derivative with respect to tjj = t y —t x . 

As for the Lorentz transformations, it is well known that an infinitesimal 
Lorentz transformation acts on a four-vector and on a spinor ip in the fol- 
lowing way: 



X» + t% x v 



(44) 
(45) 



where + e Uil = 0. By means of the fundamental relation (|3T|), we extend (|4 
to our new representation, obtaining: 



ip'^ 



o 1 



7 'W7 + (£>7 7 i y i ) 



(46) 



We then separate the boosts from the spatial rotations rewriting (46) as follows 

p = ip + ie R • Kip + i e B • B ip (47) 
where the vectors eR, and ee are 



£23 
£31 
£12 



eoi 
£02 
£03 



(48) 



and the components of the two vector operators R and B are expressed by 



R k ip 



a k 
a k 



a k 

CTfc 



(49) 



B kP = -^oi k ip - -ipa k 



(50) 
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with afc given by (p5j). 

Now we are able to find solutions of the motion equations representing particles 
moving with arbitrary speed v = vn. Let's call ifQ one of the four solutions ([36])- 
(|39|) representing particles at rest. The first step is to solve equations (|4^) and 
( f43| ) with respect to tp, obtaining a complete solution ^(xs^d) independent 
from xd and t$; the solutions (|36|)-(|39|) are recovered in the case tn = 0. Then 
we apply a boost with speed v along the unit vector n, which corresponds to the 
following transformations: 



t' s = ts cosh £ + (x§ • n) sinh £ 

x' s = x§ + (xs • n) n (cosh £ — 1) + n ts sinh £ 



(51) 
(52) 
(53) 



where the parameter £ is defined by tanh£ = v; obviously to and xq transform 
like ts and xs- Thus we obtain the following solution: 



^ v (4*Siti?,*D) = e l?n - B <A)(xs,*D) 



(54) 



It is clear that in the "physical" case t' D = x^ = the solution (54) vanishes 
outside the region defined by x s = vi^.; therefore it represents a particle moving 
with speed v perfectly localized in the position Xg=vi' s . 

Now we will express in our new representation the physical quantities associ- 
ated to the Dirac particles, which are energy, momentum, charge and spin. Let's 
start from the charge Q, whose space density in the original representation was 
the time component of the four-current 



( x ) = ^ ( x ) yv?/, (x) 



(55) 



By means of the fundamental relation (|3l]) we obtain in the new representation 
the following expression 



J»{x) = Tr{^(x s ,x D ) 7 V} 



Xs=X , Xjy=0 



(56) 



The time conservation of the charge Q is a consequence of the continuity equation 
d^J^ = 0, which is true in both representations. The explicit form of the charge 
Q in the new representation is simply given by 



f 4 
Q = / d 3 x fkk (xs, x D ) 

J fc=i 



x s =x,x D =0 



(57) 



In the same way, starting from the energy-momentum tensor in the original 
representation, we easily obtain the following expression for the same tensor in 
our new representation: 



T% (x) = Tr | -i-^- ip (x s , x D ) 7 V 



(58) 



Xa=X , Xr>=0 
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In both representations the conservation of energy and momentum is expressed by 
the continuity equation d^T^ v = 0, while the definitions of energy and momentum 
in the new representation are 



E 



S(x) d 3 x 



p % = m( x )d 3 x 



with the spatial densities £ and V% given by 



£ 



Vi 



i oo 



i0i 



Tr < i a k 



dip 
dx 



+ m (3 ip 



E 

k=l 



dx 



D,i 



D,k 



x s =x , x D =0 



x s =x , x D =0 



(59) 



(60) 



(61) 



Finally, from the angular momentum tensor in the original representation, we 
obtain the following expression for the same tensor in the new representation: 



M» Xv (x) = x A T» v (x)-x v T 



+ Tr |^99(x s ,xd)7°7 



(62) 



In both representations the angular momentum conservation is expressed by the 
continuity equation d^M ^ v = 0. Besides the angular momentum spatial densities 
are given by: 



M 0ij (o 



^ijk 



(C k {x) +S k (x) 



h ji k — 1) 2, 3 



(63) 



where e^-fe is the rank three antisymmetric tensor; C is the orbital angular mo- 
mentum density, while S is the spin density. The spin components are then given 
by: 



S t = J Si(x) d 3 x = |d 3 xTr| T 
From expressions (|56|), ( |58| ) and (|6^), 



o. L 
a; 



V?(x s ,x D ) 



(64) 



/ x s =x,x D =U 

it is clear that all observable quantities 
depend on the value of the field ip (and its space-time derivatives) in the region 
where xd = 0, or x = y. Once again we are induced to consider xs as the 
"physical" space-time coordinate, while xd is an "auxiliary" coordinate, which 
has observable effects only around the point xd = 0. 

We are now able to associate to the solutions defined in (|36|)-(|39|) their physical 
properties. From the above definitions, we see that ipA, Pb, ¥c and ipp have 
energy E = m, momentum P = and spin components S\ = S2 = 0. As for the 
charge Q and the third spin component S3, we have the following situation: <pa 
has Q = +l and S3 = +1/2, thus representing a spin up positron; ips has Q = +l 
and S3 = —1/2, thus representing a spin down positron; ipc has Q = —1 and 
S3 = —1/2, thus representing a spin down electron; finally, ipD has Q = — 1 and 
S3 = +1/2, thus representing a spin up electron. Besides, if we apply a rotation or 
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a boost to these solutions, we are certain that the physical quantities associated 
to the new solutions will have the usual relativistic values; this is true because we 
have expressed all the observable quantities in a covariant form. 

Until now, we have shown that our new representation admits solutions which 
can be interpreted as relativistic point-like particles moving at constant speed in 
an arbitrary direction, with the correct values for energy and momentum. In their 
rest frame, the particles have spin S = 1/2, while the spin orientation is arbitrary; 
as for the charge Q, we have solutions with Q = +1, to be interpreted as positrons, 
and solutions with Q = — 1, to be interpreted as electrons. By adding together 
an arbitrary number of such solutions, we obtain a system of non-interacting 
particles. Although attractive, this description is far from complete: indeed, there 
are a lot of solutions for which such particle interpretation is not possible and our 
model does not explain why the physical solutions should be acceptable while the 
non-physical solutions should be discarded. 

At this point, we must put forward the following consideration: real particles 
interact with each other and we can observe them only because they interact 
with our measuring devices; therefore, all mathematical models describing non- 
interacting particles are equally useless for explaining physical phenomena. To 
complete our model, then, we need an interaction between particles; we will choose 
the electromagnetic interaction, which has provided the first striking application 
of quantum field theory, as well as the first example of a modern gauge theory. 

To introduce the electromagnetic interaction in our model, we return to the 
original representation, where we make the usual substitution 

tt— "o ieA u (65) 

If we insert ( |65[ ) in the non-covariant Dirac equation (^0|) , and then switch to our 
new representation by means of ([H|), we obtain the following time evolution: 

= H (/? + eHi (66) 



where Ho is still defined as in (p2|), while Hi is given by 

Hi = A°(y) ip - A°{x) <p + A fc (x) a k ^ - A k (y)<pa k (67) 

where, as usual, the implicit sum is over k = 1,2,3. 

As for the covariant equations ( |40| ) and (f4l|), in the presence of an external 
electromagnetic field they become: 

i^(J^-ieA^x)^-m V = U (68) 

* (^ +ie A^))^VV + m997 = (69) 
As in the free field case, fl69^) can be obtained by taking the adjoint of (38) and 



by exchanging x with y. The non-covariant time evolution ( p6| ) can be recovered 
from (^|) and ( |69| ) simply by setting t x = t y = t. 
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Let's spend a few words about charge conjugation and gauge invariance. In 
the original representation, if ip is a solution of the Dirac equation interacting with 
an electromagnetic field 

{df.-ieA^ip -rrnp = (70) 

then the field 

tf[x) = j 2 ip*(x) (71) 



is a solution of the equation obtained from (|70|) by changing the sign of the inter- 
action parameter e: 

i-fidp + ieAJi// -mil/ = (72) 

The transformation ( [71] ) is called charge conjugation and depends on the particular 
representation ( |35| ) we have chosen for the Dirac matrices. It is easily seen that ip 
and ip' have the same charge Q = J ip^ (x) ip (x) d 3 x ; we would expect the charge 
conjugation to change the sign of the charge, as it happens for instance in the 
case of the Klein-Gordon equation, but this is of course impossible for the Dirac 
equation, since the charge is always positive in the original representation. If 
we now switch to our new representation, we find that the transformation fl7l| ) 
becomes 

y) = --/ 2 ip T {y, x) 7 2 (73) 

Again we find that ip' satisfies the motion equations obtained by changing the 
sign of the interaction parameter e, and again ip and <p' have the same charge, 
i.e. Q' = Q. However, if we now change the sign of (|73|), we still have a solution 
of the same motion equations, but now we have Q' = —Q. Therefore in our new 
representation the correct definition for the charge conjugation is 

ip'(x,y) = + 7 y r (2/,z)7 2 (74) 

It can be shown that (|74| ) changes the sign of the whole four-current J^ 1 , so that 
the field if/ satisfies the conjugate Maxwell equations for the electromagnetic field; 
this property was obviously not true in the original representation. Besides, the 



transformation (74) leaves unchanged all other observable quantities, namely en- 
ergy, momentum and spin, and this is exactly what we expect from a well defined 
charge conjugation transformation; for instance, if we consider the four solutions 
(|36|)-(j39|), we easily see that ipu (<pc) is obtained by charge conjugation from ipA 
(<Pb) and viceversa. 

Turning now to the gauge invariance, we remember that in the original repre- 
sentation a gauge transformation is given by a local phase rotation on the field ip 

iP'(x) = e id &iP{x) (75) 

together with the transformation 

^{x) = A,{x) + \^9{x) (76) 
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for the electromagnetic field. Gauge invariance means that if and are so- 
lutions of the motion equations, the same is true for t// and A'^. In our new 
representation (75) becomes 



<p'{x,y) =e i ^- e ^<p(x,y) (77) 



and it is easily seen that the motion equations (pq), or (J69|), are invariant with 



respect to the transformations (76) and (|77j), i.e. they are gauge invariant. 

Until now we have considered the electromagnetic field as an external source, 
rather than a dynamical component of our model. If we want to include the field 
A^ in our model, we have to define its time evolution in interaction with the 
"material" field cp. The dynamical behaviour of the electromagnetic field will be 
described by the classical Maxwell equations: 

d» (d^A u - dvAp) = dPFp, = J v (78) 

We simply have to define the current J v as a function of the field (p. To obtain 
the electric current density from ( |56| ) we multiply by the unit electric charge e, 
obtaining: 

J"(x) = eTr{^(x s ,x D ) 7 V} (79) 

The Maxwell equations (|78|) are again gauge invariant: the left-hand side de- 
pends only on the gauge invariant tensor F^ u , while the right-hand side depends 
only on the values of <p for xd = 0, or x = y, which are not affected by the 
transformation (|77|). 

In conclusion, we have built a model which describes the "material" Dirac field 
ip together with the electromagnetic field A^. The model is defined by the time 



evolution of the field ip, (|3g) or (|69|), and by the Maxwell equations (78) where 
the current is given by ([79]). The field ip depends on two space-time coordinates 
xs and xd', however, we point out that the duplication of the time coordinate 
is needed only for formal requirements: the motion equations (^) or ( |66| ) can 
be solved with respect to the "physical" time (ts = t and to = 0) without even 
mentioning a second time coordinate. On the contrary, the electromagnetic field 
A^ depends on a single space-time coordinate x and its time evolution is given by 
the classical Maxwell equations; it is interesting here to note that some authors 
(see Marshall and Santos |3J]) reject the quantum concept of an electromagnetic 
field made up by point-like massless "photons" , and believe that all experimental 
data involving light fields can be explained by means of the unquantized Maxwell 
equations. 



4 Discussion 

Let me explain the basic ideas (or maybe prejudices) on which the present paper 
is built. The main idea is that quantum mechanics cannot be a fundamental 
theory: I cannot force myself to believe that probabilities are an objective feature 
of the physical world, and therefore I maintain the old-fashioned belief that the 



14 



fundamental laws of nature must be deterministic and quantum mechanics must 
be just a statistical formulation of these fundamental laws. The obvious example 
for clarifying this concept is the classical motion of a free non-relativistic particle, 
with position x and momentum p : the "wave equation" 

%=--¥■ < 8o > 

at m ox 

for the joint probability distribution / (x,p) is just a statistical formulation of the 
deterministic state equations 

x = — p = (81) 

m 

We are not worried about the "collapse of the wave function" at the time when 
the position x is measured, because we do not interpret / (x,p) as a real objective 
field; in any case, such interpretation would be strongly hampered by the fact that 
/ (x,p) is defined in configuration space, whose dimension depends on the number 
of particles described. 

Now, what mathematical form should have this deterministic law which will 
finally replace quantum mechanics and quantum field theory? I am convinced 
that it should be a classical (unquantized) field equation, satisfying the principle 
of relativistic covariance. However, since until now all efforts at finding this uni- 
fied field equation have failed, we probably need to extend the classical concept 
of field, to gain some room for maneuver. In my opinion, the most natural ex- 
tension is to allow for extra space-time dimensions, following the idea which was 
originally proposed by Kaluza and Klein [^], and has been recently revived by 
Hasselmann ||. Another possible choice could be to abandon the concept of a 
continuous space-time, and suppose that the ultimate structure of space-time is 
discrete (see for instance Budnik ||); however, besides having a prejudice against 
"objective probabilities", I also have a prejudice in favour of Newton's "Natura 
non facit saltus", basing this prejudice mainly on the wish to maintain the beau- 
tiful space-time simmetries and conservation laws which are strictly connected to 
the continuous structure of space-time. 

Thus, we are looking for a relativistic field equation (or a system of equations) 
defined on a space-time with more than four dimensions. But what should be 
our starting point for this long and perilous search? The authors cited above 
stated immediately their ambitious goal: they wanted to unify gravity with all 
other known interactions (in the Kaluza-Klein case it was just the electromagnetic 
interaction, while Hasselmann's "metron" model includes all Standard Model's 
interactions). Even if this is indeed the final goal, I propose here a more cautious 
approach: instead of trying to guess from the beginning the final form of our field 
equations, let's start from the field equations on which quantum mechanics and 
quantum field theory are based (the non-relativistic Schrodinger equation and the 
relativistic Dirac equation) and let's try to extend them in such a way that their 
solutions may now be interpreted as real objective fields. 

In the introduction to this paper (Section ||) I listed the two main obstacles to a 
realistic interpretation of the wave function in non-relativistic quantum mechanics: 
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(1) the wave function spreads out with time and (2) it is defined on a variable 
dimension configuration space. Well, in Section ^ we saw that if we choose the 
Schrodinger equation for density matrices, instead of the usual equation for pure 
states, it is possible to find solutions which are perfectly localized both in position 
space and in momentum space. In the free field case, these solutions behave 
exactly as non-relativistic point-like particles moving at constant speed, with the 
correct values for the observable quantities (position, energy and momentum). By 
adding together an arbitrary number of such solutions we obtain a field describing 
a system of non-interacting particles; this field depends always on the same two 
space coordinates, the "physical" coordinate xs and the "auxiliary" coordinate 
xr>. Both problems (1) and (2) are solved, and I believe that this is a first (small) 
step in the right direction. 

It is interesting here to note that some authors already consider the Schrodinger 
equation for density matrices more fundamental than the one for pure states. 
For instance Olavo || explicitly names the equation for density matrices "First 
Schrodinger's Equation", while the equation for pure states, named "Second 
Schrodinger's Equation", is obtained from the first in the less general case where 
the density matrix is a product of pure states. Olavo also defines the position 
operator and the momentum operator for density matrices and shows that they 
commute, similarly to what we did in Section ^; from this property Olavo infers 
the "negation of ontological origin of Heisenberg's uncertainty relations" and I do 
indeed agree with him. However, Olavo's purpose is to demonstrate that quan- 
tum mechanics can be derived from classical mechanics, and therefore he is not 
interested in giving a realistic interpretation of the density matrix. 

A second example can be found in Mermin's "Ithaca interpretation" (7j; in this 
paper the author lists six desiderata for a satisfactory interpretation of quantum 
mechanics, and then writes "If you take Desideratum (5) seriously (i.e. generalized 
Einstein locality), then there can be no more objective reality to the possible 
different realizations of a density matrix, than there is to the different possible 
ways of expanding a pure state in terms of different complete orthonormal sets 
... In the case of an individual system the density matrix must be a fundamental 
and irreducible objective property, whether or not it is a pure state". This last 
statement sounds quite similar to the idea that I am suggesting in the present 
paper, i.e. the density matrix as a real objective field. However, this resemblance 
disappears when we read Desideratum (6), which states that probabilities are 
objective intrinsic properties of individual physical systems; unlike Mermin, I still 
believe that probability is "just a way of dealing sistematically with our own 
ignorance" , and not a fundamental feature of the physical world. 

Let's turn now to the relativistic Dirac equation. We immediately see a new 
obstacle to a realistic interpretation of its solutions: (3) the charge density ip^ip is 
always positive, and therefore the classical Dirac field does not allow us to describe 
particles with both positive and negative charge, i.e. positrons and electrons. The 
second quantization procedure solves point (3) by means of the anti-commuting 
properties of the creation-annihilation operators, but the points (1) and (2) listed 
above remain unsolved; furthermore, the final wave function (or state vector) is 
defined in a Hilbert space whose dimension is dramatically increased. On the 
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contrary, in Section |3] we saw that, by extending to the Dirac equation the same 
representation already used for the non-relativistic Schrodinger equation, it is 
possible to find solutions which behave as relativistic point-like particles moving 
at constant speed, with the correct values for all observable quantities. By adding 
together an arbitrary number of such solutions we can describe a system of non- 
interacting particles, and the resulting field depends always on the same six space 
coordinates, three "physical" coordinates xs and three "auxiliary" coordinates xd- 
By introducing also two time coordinates (a "physical" time coordinate ts and an 
"auxiliary" time coordinate t^) we can recover the formal Lorentz covariance of the 
Dirac equation. Finally, in the new representation the charge can have both signs, 
and more generally it is possible to define a charge conjugation transformation 
which changes the sign of the charge and leaves unchanged all other physical 
quantities. Therefore, points (1), (2) and (3) are all solved, and I think that this 
is a second step in the right direction. 

Thus, for both the non-relativistic Schrodinger equation and the relativistic 
Dirac equation we were able to find solutions which can be interpreted as non- 
interacting point-like particles. At this point, however, we must face a big problem: 
these solutions are only a limited subset of all possible solutions, while for most 
solutions this particle interpretation is not possible; until now our model does not 
explain why the "particle" solutions should be accepted while the other solutions 
should be discarded. I will discuss this problem only for the Dirac equation, ne- 
glecting the non-relativistic case for the following reasons: first, not all potentials 
V (x) have physical meaning; second, since non-relativistic quantum mechanics is 
just an approximation for relativistic quantum field theory, if we find an accept- 
able model for quantum field theory then the corresponding model for quantum 
mechanics will follow in the non-relativistic limit. 

In my opinion, there is only one way to single out the "particle" solutions of 
the Dirac equation from all other solutions: we must add a non-linear interaction 
to the model. The interaction must be non-linear for obvious reasons: for instance, 
the field obtained by changing the sign to a positive energy solution cannot be 
itself a solution, because it would have negative energy (negative mass in its rest 
frame) and therefore it would be non-physical. Thus, in the second part of Sec- 
tion |H we inserted in our model the electromagnetic interaction, by means of the 
usual "minimal coupling" prescription. We also included in our model the dy- 
namical evolution of the electromagnetic field, in the form of the classical Maxwell 
equations; the four-current in the Maxwell equations depends on the Dirac field 
through equation (79). Our final model provides a complete description of the 
mutual interactions between the "material" Dirac field and the "gauge" field, i.e. 
the classical electromagnetic field. 

Now the question is: does this final model really solve the problem of elimi- 
nating the non-physical solutions? This will be true if the only stable solutions 
for the "material" field tp are the vacuum (i.e. the solution (f = 0) and the single 
particle solutions, which must be strongly localized in physical space and must 
provide the correct values for all observable quantities; we can relax this last re- 
quirement by accepting some form of renormalization, i.e. the physical mass and 
charge m p h and e p h could be different from the "bare" quantities m and e which 
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appear in the model equations. Unfortunately, the task of finding the solutions 
to our model (either by symbolic calculations or by computer simulations) seems 
to me extraordinarily complex, so our question will find no answer in the present 
paper. However, let's suppose that the answer may be affirmative; if this were 
true, then we could attempt to interpret our model as a deterministic alternative 
to quantum electrodynamics. To reach this ambitious goal, we would have to sat- 
isfy (at least) two further requirements: first, in scattering processes our model 
should be able to reproduce the well known phenomenology of QED, with the 
same probability amplitudes; second, by adding an external source to our model 
(for instance the charge of the hydrogen atom nucleus), the stable solutions should 
reproduce the same energy spectra which can be obtained from standard QED. 

On the contrary, if the answer to our question were a plain no, this would not 
necessarily entail that all our previous efforts have been useless; we could try to 
make some changes in our model, along one of the following lines: 

• Maybe the electromagnetic interaction alone is not enough to single out the 
"particle" solutions from all other solutions: maybe it is necessary to insert 
in our model the other known interactions, or even some new (yet unknown) 
interaction, which would be responsible for the strong localization of the 
electron and positron fields. My guess that the electromagnetic interaction 
may suffice is based only on the fact that QED can exist as an independent 
theory even if we know that it is really a part of a more general theory, i.e. 
the Weinberg-Salam model for electroweak interactions. 

• In this paper we suppose that the "auxiliary" space-time, indexed by the 
coordinate xd, is infinite and flat. It could as well be finite and strongly 
curved, similarly to the fifth dimension in the Klein version of the original 
Kaluza theory. If the curvature radius were really small compared to macro- 
scopic space-time distances, this could explain why we do not observe the 
extra space-time dimensions. 

• Finally, we could look for a way to extend the electromagnetic field to the 
extra space-time dimensions, obtaining a new field (xs,xd) which would 
be equal to the classical electromagnetic field in the case xd = 0. The 
resulting model should have a better simmetry between the material field 
and the gauge field. 

Of course, I have no absolute certainty that the method presented in this paper 
will finally produce acceptable explanations of physical phenomena: I just tried to 
find a new approach to deal with a very old problem. However, it can be shown 
that if we apply this method to simple non-physical systems (two- and three- 
dimensional Hilbert spaces) we indeed obtain deterministic models equivalent to 
the corresponding quantum systems; in this case, the density matrix representation 
becomes simply the adjoint representation of the SU(2) and SU(3) groups. 

In conclusion, I will point out some features of our model which could be re- 
lated to the debate about quantum non-locality; let's briefly summarize the main 
results of this long lasting debate. In 1964 Bell |8| showed that no local hidden 
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variable (LHV) model can reproduce the same probability correlations obtained 
from quantum mechanics; he did so by constructing an explicit relation (the first 
famous Bell inequality) which must be satisfied by every LHV model and which 
is violated by some quantum states. In the early 80's, experiments were per- 
formed |9], [n]] which seemed to indicate that Bell inequalities are indeed violated 
by nature. If this were true, then we would be forced to choose between determin- 
ism and locality; for instance, Mermin Jjj] chooses to abandon determinism and 
to retain locality. However, the debate about the experimental violation of Bell 
inequalities is still open: it seems that this violation can be inferred from the exist- 
ing experimental data only by introducing some supplementary assumptions (no 
enhancement, faithful sampling, . . . ) and therefore many authors [11, [l^] believe 



that until now no incontrovertible violation has been really observed. Besides, 
Marshall and Santos [||, 13] are convinced that the corpuscular interpretation of 



light plays an important role in the above-mentioned assumptions, and were able 
to build a local "wave" theory of light, derived from Maxwell's classical equations, 



whose predictions are in good agreement with the experimental data [14]. 

The model presented in Section [3| of the present paper is expressed in a rel- 
ativistic covariant form and therefore it must satisfy the principle of relativistic 
causality (space-like separated events cannot influence each other and no interac- 
tion can propagate at superluminal speed). However, the following features seem 
to indicate some form of non-local behaviour: 

• Expression (^) establishes an instant relation between the fields <p (x, y) 
and </?(y,x). However, this non-local relation is confined to extra-space, 
while in the physical case x = y we obtain a local hermiticity property for 
the matrix field ip (x,x). 



The solutions (|36|)-(|39D representing particles at rest are perfectly localized 
in the physical position xg = and do not depend on the auxiliary coordi- 
nate xd- However, it can be shown that the solutions representing moving 
particles, while being perfectly localized at xs = vt when xd = 0, spread 
out along the x§ axis when xd 0. Indeed, if n is the direction of motion, 
then the corresponding field is different from zero inside the region defined 
by [xg-n — vt\<v |xd-n|/2. Therefore, if we have two such solutions local- 
ized at two distant physical points, their fields will anyway overlap in some 
extra-space regions. This does not mean that the two fields interact, since 
in the free field case the motion equations are linear, but in any case this is 
again a form of non-locality in extra-space. 



Finally, equation (67) clearly shows that the electromagnetic field M 1 (z) acts 
instantly on the field ip (x,y) at all those points having x = z or y = z. There- 
fore an interaction can propagate between two physical points A = (za, za) 
and B = (zb,zb) through the following mechanism: the field ^(za,za) 
acts instantly on the electromagnetic field (za) through the current (|79|); 
the field A^ (za) acts, again instantly, on a extra-space point C = (za, zb); 
finally, the interaction propagates from C in the x direction until it reaches 
the point B = (zb,zb), this time moving at the speed of light. The time 
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elapsed is exactly the time needed for a light signal to travel from A to 
B, but the interaction has been propagated entirely through extra-space: 
at intermediate times, the effect of the interaction is not observable at any 
physical point between A and B. This seems to be a non-local behaviour 
with physical effects, and it would be interesting to investigate if this kind of 
non-locality is in some way related to the non-local features emerging from 
quantum mechanics. 
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